Ultrasensitive Atomic clock with single-mode number-squeezing 



L. Pezze^ and A. Smerzi^ 

^ Laboratoire Charles Fabry de I'Institut d'Optique, 
campus Polytechnique RD128, F-91127 Palaiseau cedex, France 
^ CNR-INFM BEC center and Dipartimento di Fisica, Universita di Trento, 1-38050 Povo, Italy 

(Dated: May 3, 2010) 

We show that the sensitivity of an atomic clock can be enhanced below the shot-noise level by 
initially squeezing, and then measuring in output, the population of a single atomic level. This can 
simplify current experimental protocols which requires squeezing of the relative number of particles 
of the two populated states. We finally study, as a specific application, the clock sensitivity obtained 
with a single mode quantum non-demolition measurement. 
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Introduction. Atomic clocks and interferometers are 
among the most sensitive measurement devices avail- 
able within the current technology [Ij. Their precision 
is bounded by the fundamental noise imposed by quan- 
tum mechanics uncertainties. With uncorrelated atoms, 
the phase sensitivity of the Ramsey sequence scales as 
the inverse square root of the total number of particles, 
the so-called shot noise (or standard quantum) limit, first 
observed experimentally in Q. The possibility to over- 
come the shot noise limit by quantum engineering spe- 
cific atomic correlations is a break grounding prediction 
which is under intense experimental investigation. Most 
current efforts focus on the creation of spin-squeezed 
states 0-01 ■ This can be achieved, for instance, by ma- 
nipulating a cloud of cold atoms via the back-reaction 
of quantum non-demolition (QND) measurements 0, [E| 
or with interaction-induced nonlinearity using Bose Ein- 
stein Condensates (BECs) With BEC, entangled- 
enhanced Ramsey phase sensitivity has been recently 
demonstrated [1]. 

In the current literature, it is shown that sub shot noise 
phase sensitivity in an atomic clock is generally associ- 
ated to spin squeezing. In this manuscript we demon- 
strate a sub shot noise phase sensitivity up to the Heisen- 
berg limit with the initial squeezing and the measurement 
in output of the particle population of a single atomic 
level. This can simplify current interferometric proto- 
cols. A non-destructive atom-light interaction of a sin- 
gle clock level has been experimentally demonstrated in 
[1, 1^. Therefore, our prediction can be readily tested 
experimentally and find application for precision atomic 
sensors within the present state-of-the-art technology. 

Sub shot-noise with single mode squeezing. We consider 
two atomic clock levels a and b of energy hwa and fiwb, 
respectively. The goal is to estimate the frequency dif- 
ference Aw — uja — ujf, with the highest possible sensitiv- 
ity. The Ramsey interferometric sequence consists of four 
steps: a Rabi 7r/2 pulse of constant power and frequency 
a; (as close as possible in resonance with the atomic tran- 
sition) applied for a time r = 7r/2r2ij, being the Rabi 
frequency. Then the system freely evolves for a period 
of time T. Finally, after a second 7r/2 pulse, the num- 
ber of particles is measured in a single output level. The 
quantum mechanical expectation value of the number of 
particles in the output a mode is given by ( 

^a)out — 



(na)inpCOs2(f)-f(n6)inpSin^(|)-i(a'f'&+S'''a)i„pSin6'[i3l, 
where a {a)) and b [w) are particle annihilation (cre- 
ation) operators of the a and h mode, respectively, and 
fia = a^a {fib = b^b) is the number of particles operator. 
The quantity {fia) out depends on the phase 9 = 6 xT ac- 
cumulated during the free precession, being S = Auj — uj 
the detuning of the Rabi pulse from the atomic transition. 
By collecting m measurements with results nl^\ 
we can calculate the average number of particles in the 
a output mode n°"* = "-o V™- The phase is in- 

ferred by approximating the expectation value (a^a)out 
with the classical average and inverting the equation 



(a^a) 



,cos2(%t) + (5t5).„pSin2(%t) 



&^a)inpSin6'ost, where ^est is the estimated value of 6*. In 
the central limit, the expected phase sensitivity is 



A9 = 



(Ana)c 



m |d(nQ)out/d6'| 



(1) 



where A9 is the variance of ^cst, calculated with error 
propagation, assuming that fluctuates with variance 
(Ana)out = ("-a)out " (?^a)ouf In particular, when 6* 
and with initial symmetric populations (na)inp — (^fc)inp, 
Eq.([T]) becomes 



A9 



2{Ana)i 



TO|(at&-|-6ta);„p| 



(2) 



Equations ([T]) and © relate the interferometric phase 
sensitivity A9 to the fluctuations of the number of par- 
ticles in a single input atomic level, and capture the 
main results of this manuscript. In particular, accord- 
ing to Eq. ^ , sub shot- noise phase sensitivity can be ob- 
tained by squeezing the population fluctuations in the 
input a, {Ana)fnp < while keeping the coher- 

ence between the two mode (a^b) 



inp 



(n)/2. Here 

{n) = {riajinp + {n-b) inp is the average number of par- 
ticles in the input state. As an example, we con- 
sider the product of two Gaussian pure states \ip)a,b 
I]^^exp[-(n- (n)/2)2/4cr^_J|n)a,h as input of the in- 
terferometer, where ab = ^{h)/2 and CTq = i^y^ {n) /2. 
Number-squeezing in the a mode [TTj is obtained for 
K < 1, and equation ^ predicts a sub shot noise sensi- 
tivity A9 = K/y/m{h). In the Fock limit (7a ^- 0, Eq.© 



2 



predicts the ultimate Heisenberg scaling 



(3) 



The Fock state limit is particularly interesting and will 
be further discussed below. 

Notice that Eq. ^ resembles the familiar relation be- 
tween interferometric sensitivity and spin squeezing jl2l - 
[T3 | . However, since the total number of particles is not 
fixed, number squeezing in a single clock level does not 
necessarily implies the squeezing of the relative popula- 
tion between the two input modes. Only when the total 
number of particles is fixed, the two become equivalent. 
Our analysis includes this special case. 

Entanglement and Fock state limit. The relation 
between multiparticle-entaglement and sub shot noise 
phase sensitivity in a linear interferometer has been re- 
cently discussed for states of fixed [iHl and fluctuat- 



17| number of particles (qubits). In particular. 



mg 

states satisfying the inequality Fq [p, Jy] > (n) , being 
Fq [p, Jy] the quantum Fisher information jl^, and 
Jy = {a)h + Wa)/2 [iGj, are entangled and provide a sub 
shot-noise phase sensitivity in a Ramsey interferometer. 
The optimal phase sensitivity, after m independent mea- 



surements, is lS.6 = l/y mFQ[p, Jy\. In particular, for an 
input state 

Pinp=\N)a{N\® pb, (4) 

being p^ — P'".W)b{n\ a generic density matrix of 
= Yl,n=o P^n average number of particles, wc find 



M = 



1 



my/2N (hb) inp + N + {hb) inp 



(5) 



Equation ([5]) is characterized by interesting limits. If 
the mode b is left empty, pn = Sn,Oy the sensitivity of 
the clock is given by the standard quantum limit (SQL), 



Ae 



where 



N 



[rib 



is the aver- 



age number of particles in the input state (in this case 
(fi) = N). Conversely, if the mode b is left almost (but 
not completely) empty, with (nb)inp ^ the sensitiv- 
ity A9 « l/y/2{hb)inp 



a factor ^2(nb)i,ip 
dition N = {fib) 



ly/m{h) is below the SQL by 
1. Moreover, at the optimal con- 
Eq.® predicts AO ^ ^ 



dependently of the input states pb, and, in particular, 
we recover Eq.Q. Measuring the number of particles in 
a single output port not only provides an optimal esti- 
mation strategy but also saturates the Heisenberg limit 
17| . As anticipated above, we note that Eq.Q is not 
a spin squeezed state. The spin squeezing parameter 
e = {AJ.)lp/{{J.)lp + {Jy)lp) Mm diverges in gen- 
eral when calculated for the input state Eq.(|3]), and is un- 
determined, = 0/0 , when also the input 6 is in a Fock 
state {pn = Sn,N) UM- the latter case, corresponding 
to the twin Fock state, Eq.® recovers the Heisenberg 
limit sensitivity first predicted in [20l |. 

QND state preparation. We now discuss a possible 
experimental implementation of the protocol discussed 



above. The single-mode number squeezing is produced 
by a QND interaction between the atoms cloud and a 
light field, recently demonstrated in [sl, 0. The ini- 
tial atomic cloud, obtained by optically pumping the 
atoms to a single level, is described by the density ma- 
trix p = En%Pn \N,0){N,0\, where |7V,0) = |iV),|0)fc 
is the state with N atoms in the levels a and atoms in 
level 6, which occurs with probability P/v (Pn > and 
X)n Pn = !)■ In order to equally populate in average the 
two input modes, we apply a 7r/2 pulse. The density ma- 

Y,'n^qPn\^n){iPn\ where \i>N) = 
n\{N-ny. 1'^'^ ^ The number of parti- 
cles distribution in each mode has (ha) — (fib) — {n)/2, 
{Ahaf = {Aubf = (a2 + (n))/4 where (n) = Y^nPn^ 
and cP' = X^AT Pn{N — {n)Y . With this state as input of 
the clock, Eq.dl]) predicts 



trix becomes p"^ 

S^N 1 



Ae 




cos^ 6 



(6) 



At the optimal value of the phase shift 6 — it/2, the 
phase sensitivity is at the standard quantum limit AO = 
1 / yjm{n). It is possible to overcome this limit by squeez- 
ing the number of particles fluctuations in one mode. 
This is done here by letting the atoms in a to inter- 
act with a coherent light field |a) (we take, without loss 
of generality, the amplitude a to be real) via the QND 
Hamiltonian 

i?QND = % (c^c) {a) a), (7) 

where & (c) created (annihilate) a photon of light field 
mode and g is the coupling parameter. The Hamiltonian 
Eq.© describes an out of resonance interaction of the 
light field with the atoms in the level a [ij]. We assume 
here that the detuning of the light is large enough to 
adiabatically eliminate the excited state population and 
neglect possible decoherence effects [2l[. The coupling 
strength is normally weak but can be enhanced by plac- 
ing the atoms inside an optical cavity [23 |. The entan- 
gled atom-light system is described by the density matrix 
= e^HQ^-Dt/h^m ^ |Q,)(Q,|e-»-^QNDt/n^ ^jjgj.g ^ jg ^i^g 

interaction time between the light and the atomic cloud 
and a is the amplitude of the coherent state. The effect 
of the QND interaction is to shift the phase of the light by 
a quantity proportional to the number of atoms. There 
are different experimental possibilities to estimate the 
phase shift [I,!!]. In the following we consider the homo- 
dyne measurement of the p quadrature, p = {c— ct)/2i. 
The probability to measure the eigenvalue p is given by 



^(0)1 



is traced 



Poip) = Tr[p^"V)(p|J, where p^"^ = TiaIp^alI 
over the atomic degrees of freedom [2^. After the mea- 
surement, with result po, the density matrix of the atomic 



cloud becomes 



"(1) rp 

Pa = Tri 



Pal bo) (Po I 



/Po(po). The 



back reaction effect induced by the measurement [2J| re- 
duces the fluctuations of the number of atoms in the a 
mode. In the limit n{h) <^ 1, with D, — gt, the quadra- 
ture probability is a sum of Gaussian functions [25| and 
we can distinguish two regimes. When aVL ^ 1, the mea- 
surement result is only compatible with a precise number 
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of atoms and therefore the QND interaction projects the 
atomic population of the a mode to a Fock state. In 
this case the sensitivity is at the Heisenberg hmit Eq. (jS]) . 
Conversely, when aO. < 1, which is the most realistic 
scenario from the experimental point of view, the result 
of the quadrature measurement is compatible with sev- 
eral values of the number of particles in the a mode, and 
the back reaction only produces a moderate squeezing. 
Therefore, in order to increase the squeezing, it is neces- 
sary to repeat M ^ 1 times the QND protocol [l^ . Ana- 
lytical results can be obtained in the limit 7 = a^^^M < 

1. In this case we obtain (Aa'^a)^ = }" t^l^]y~\ and 

{a)h + &^a)inp ~ 2(n) cos(Ma^r2). Conversely, the input 
6 retains, in average, the initial number fluctuations. By 
tuning Ma^g to values close to an integer multiple of tt, 
Eq.© gives 



A9 = 



1 



{h)y/TnV 1 + '-f{a^ + (fi)) 



(8) 



Equation (|5]) is obtained for 9^0 and arbitrary fluc- 
tuations a. In particular, when 7 > a'^/n{n + cr^), we 
obtain a sub shot noise {A9 < 1/ yinym) limit. In figure 
([T|) we plot the phase sensitivity, calculated using Eq.([T]) 
with a coherent state having ~ (h) ~ 10^, as a func- 
tion of the phase shift 9. The different lines correspond 
to different values of 7. In particular, the dotted black 
line is obtained for 7 = and is given by Eq.®. The 
upper (lower) horizontal lines is \/ yj{n) {\/2/{h)). The 
figure shows that sub shot noise can be obtained around 
an optimal value of the phase shift, 9opt (minimum of 
each curve). In figure ([Jla) we show 9 opt as a function of 
7: it is close to Oopt ^ 7''/2 for moderate squeezing and 
rapidly tends to 9 opt ^ by increasing 7. In figure ([2Jb) 
we show the optimal phase sensitivity as a function of 7. 
Numerical results (circles) agree with the solid blue line 
given by Eq.® in the limit of relatively small 7, where 
the optimal phase shift is 9opt ~ [27]. For 7 > 1 the 
QND project the state close to Fock, and the phase sen- 
sitivity converges to Eq.® (lower solid horizontal line). 
Notice that Eq.® applies also when the total number 
of particles is fixed (a = 0). In this case, the squeezing 
of the number of particles of a single level is equivalent 
to relative number squeezing. In figure ®B) the dotted 
red line corresponding to the case cr = (and (n) = 10^ 
atoms), superposes to the dots, obtained from numerical 
simulations with = {fi). This clearly shows that rel- 
ative number squeezing does not provide any advantage 
with respect to single mode number squeezing: for the 
same number of particles and squeezing parameters we 
obtain the same level of sub shot noise interferometric 
sensitivity. 

Discussion. As discussed above, a QND measurement 
on a single energy level is sufficient to prepare an input 
states useful to reach a sub shot noise phase sensitiv- 
ity. This can potentially simplify current experimental 
schemes. The reduction of the fluctuations of the rel- 
ative population (while preserving the coherence) with 
a QND atom-light interaction has been recently demon- 
strated. This has been done by i) addressing the two 




FIG. 1: Phase sensitivity Eq.([T| as function oiO and for differ- 
ent values of 7: for -y/n = (dotted black line) the sensitivity 
is given by Eq.®, j/tt = 10~^ (dot-dot-dashed green line), 
7/7r = 10"* (dot-dashed blue line), 7/7r = 10~^ (dashed red 
line), 7/7r — 10~^ (solid black line). The curves have been 
obtained for a Gaussian distribution of the total number of 
particles (cr^ — (n)) and {h) — 10^ atoms. The upper (lower) 
horizontal lines is {n) {\/2/{fi)). 
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FIG. 2: a) Optimal phase shift 9opt (minimum of Eq.([T])) as 
function of 7. b) Phase sensitivity Eq.([T]) calculated at the 
optimal point ^optas function of 7 (dots). The lines are Eq.® 
for (T^ = (n) (blue line), and a = (dotted red line). The 



horizontal lines are 1/y 
line). 



'{h} (upper line) and Eq.® (lower 



atomic levels with two carefully detuned laser beams |j| 
or ii) with a single laser beam tuned at a very precise fre- 
quency and detuning where the index of refraction of the 
gas crosses a zero value [sj. On the other hand, when ad- 
dressing a single level, it is possible to perform QND with 
the freedom of choosing the laser frequency and detuning 
so to minimize the decoherence effect due to spontaneous 
emission. As shown in 4, 51 this is one of the main lim- 
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itations of the level of squeezing reached experimentally. 
Within our scheme, it is also possible to select the atomic 
transition in order to maximize the response of the light 
field to the atom-light interaction. A further experimen- 
tal advantage of our proposal is that the clock frequency 
can be estimated by measuring the number of particles in 
a single output port of the Ramsey interferometer. This 
avoids the further noise introduced by inverting and mea- 
suring the population of the second mode, as currently 
done. We thus expect that our protocol might exper- 
imentally provide higher squeezing with a more robust 
apparatus. 

Conclusion. We have discussed a new protocol to reach 
a sub shot noise sensitivity (up to the Heisenberg limit) 
in an atomic clock. This requires: i) reduced particles 



number fluctuations in a single input mode and ii) the 
measurement of the number of particles in a single out- 
put. Our results can be interpreted in terms of useful 
entanglement created by squeezing the population fluc- 
tuations in a single mode. We provide a simple analysis 
of the number squeezing produced by a QND interaction 
between a light field and the atomic sample showing, for 
this experimentally relevant situation, the possibility to 
readily verify our predictions. Since Fock states of a small 
number of photons are currently experimentally available 
[28l |. our predictions can be relevant also in the optical 
domine. 
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